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Abstract 

We have performed numerical simulations of dendritic growth at very low un- 
dercoolings in two spatial dimension using a phase-field model. In this regime 
of growth, the dendrites present sharp corners in the tip region while the 
trailing region is parabolic, and the corresponding side-branching structures 
resemble the shape of the tip. The scaling vp^ ~ constant., where p is the 
tip radius of curvature from the fitting to the parabolic trailing region, still 
holds approximately. We find that the values of vp^ are consistent with those 
given by microscopic solvability theory. The sharpness of the tip region of 
the dendrite can be characterized in terms of the deviation A with respect to 
an Ivantsov parabola. We observe that this length scales as A ~ p, consistent 
with experimental measurements. 
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Dendritic crystal growth constitutes one of the most interesting examples of pattern 
formation phenomena in nonequilibrium dynamical systems . While most studies 

of dendritic growth have concentrated on relatively high undercoolings, recent experimental 
measurements in dendritic growth at low undercooling have revealed a new feature of 
this nonequilibrium process, that is, the appearance of a different dendritic morphology 
with a faceted structure and sharp corners at the dendritic tip. This regime, being close to 
equilibrium, has been addressed theoretically p[ to give an explanation of the self-similarity 
of the interfacial profile for different undercooling values. There are however still many open 
questions regarding the dynamical evolution of the pattern, the shapes of the side-branches, 
the detailed study of the interfacial roughness near the dendritic tip, and, most importantly, 
whether new elements should be included in the theoretical description, or if such patterns 
selected by nature are compatible with our current theoretical understanding 

Motivated by the interesting results of Refs. and [§ , in this paper we use a phase- field 
model to investigate the small undercooling regime of dendritic growth with an emphasis 
on steady-state pattern selection. Phase-field models have proven to be a breakthrough for 
the numerical simulation of unstable interfaces [jlO|,|lT],|l2i , and so far very impressive results 



have been obtained ||TT|J13|1 . The method is especially useful here since the problem involves 
high values of interface curvature. Our numerical data show that dendrites with sharp- 
cornered tips and sharp-shaped side-branchings (see below), are the steady-state pattern at 
low driving force. With appropriate characterization of the size of the sharp-cornered tip 
(see below), certain scaling laws hold. It is also found that the relation between growth 
velocity at the dendrite tip and tip radius at large driving force, still holds here. Finally, 
our numerical data are consistent with experimental measurements M. 



We use the phase-field model proposed for solidification in a pure liquid |ri] , p!2| , which 



has been shown to provide a reliable means of studying dendritic growth phenomena [IT3 
The model includes two fields. First, a dimensionless temperature field u is defined as 
u{x,y,t) = {T{x,y,t) — Ts)/(Ts — T^), where T5 is the melting temperature of the solid 
phase and Tl is the temperature of the liquid phase far from the interface. Second, a phase 
field is the order parameter of our system; thus = and 0=1 represent solid and liquid 
phases, respectively. The interface locus is determined by positions at which = 1/2. The 



dynamical equations are the following |]13 
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where the generalized Laplacian operator C is defined as C{r}{9)) = dx{r]{9)'q' {9)dy-) -\- 
9y{v{(^)v'{(^)(^x-) + V(77(^)^V-)). The parameters A, a and m are defined as functions of 
the latent heat per unit volume L, the specific heat c, the interfacial energy and mobility 
||r3| . In particular the dimensionless undercooling A is defined as A = c{Ts — Tl)/L. The 
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interfacial width is given [0 by the length scale e. In the so-called thin interface limit 
(e 0) [0 one indeed obtains the equivalent macroscopic equations |I| from the present 
phase-field model, as shown by McFadden et al. Thus instead of solving macroscopic 
equations we shall focus on solving the phase-field model directly. 

In the above equations the function r] accounts for the symmetry of the crystalline struc- 
ture of the solid, and therefore the anisotropy of the kinetic and surface tension coefficients 
at the interface. Here, as usual we have chosen a fourfold symmetry, ri{9) = H-7cos(4^), 
where 7 is the surface tension anisotropy parameter. We have been motivated to choose this 
form of anisotropy for two reasons. First, it is the simplest form and is used in the solvabil- 
ity analysis of dendritic growth , thus comparison to that theory maybe made at least at 
large driving forces or small 7's. Second, this form is adequate for the description of the 
physics of missing orientations [lT5ip!6| , |l7| in equilibrium crystal shapes. It is well known that 
missing orientations lead to sharp corners in equilibrium crystal shapes [T^JT^, hence they 
should play a role when the growth is very slow. With this form of the anisotropy it is easy 
to obtain the equilibrium missing orientations using the Frank diagram |jl8|,[l6|. 

We have integrated Eqs. (|l|) and (||) on a square lattice of width w G [100, 200] and height 
h G [1000, 2000] using an Euler algorithm with mesh size 7.5 x 10"'^ and time step 10~^. The 



values of the parameters were a = 400, m = 0.05 and e = 5 x 10^ ||13|. We have investigated 
the evolution of the system for values of undercooling A ranging from 0.02 to 0.35. To our 
knowledge no other numerical simulations at such low undercooling have been carried out 
so far. A large computational effort is required due to the increasingly demanding time and 
space factors. For instance the computational time can be estimated to be proportional to a 
large negative power of A [^ . We have fixed the surface tension anisotropy parameter 7 = 
0.1; in equilibrium this would correspond |T^ to missing orientations of angles 6 smaller than 
±21.65°, i.e. surfaces with normal angles in this range are thermodynamically unstable and 
do not appear on the equilibrium crystal shape. Of course we are facing a non-equilibrium 
situation here, but this construction for equilibrium shapes gives us some idea of what might 
happen. An additive Gaussian white noise of intensity 10~^ has been added to Eq. (^) with 
the purpose only of helping to trigger the growth of side-branches, assuming the role usually 
played by numerical noise. However noise with such a low intensity cannot affect the tip 
profile. A detailed description of the stochastic effects on dendritic growth will be presented 
elsewhere. 

Our starting configuration was a small circular solid nucleus in the center at the bottom of 
the integration grid. We found that dynamical evolution of the unstable interface progresses 
as follows. First, there was a relatively long transient, in which the temperature evolved from 
a initial step-like profile to a diffusive one, with a diffusion length comparable to but smaller 
than the size of the system. The velocity of the interface was constantly adjusting until the 
temperature flux at the interface acquired a value compatible with the undercooling imposed. 
Then, the truly dendritic evolution began. We observed that at these very low undercoolings 
the tip was no longer parabolic, but had a polygonal structure with a sharp corner (Fig. p. 
This in some sense resembles the equilibrium crystal shape, namely a "square" with rounded 



3 



edges but still sharp corners for a fourfold anisotropy |T6[ . In our simulations we constantly 
monitored the interfacial profile to find the steady-state tip shape, which usually appeared 
at the precise moment when the first side-branch started to grow. After this moment the 
sharp-cornered tip structure changed very little with time. We note that the only important 
difference between this simulation and the previous ones [|I^ is the very small undercooling, 
and thus, very small driving force. However notice that the steady-state shape including the 
shape of side-branches is here quite different from that which appears in the high driving 
force case. Here growth of side-branches is largely inhibited as their sizes stay more or less 
constant along the sides of the dendrite. We have checked this side-branching structure by 
doubling the system size and obtained the same result. The sharp-cornered tip structure 
can be characterized using a length A which corresponds to the distance between the tip and 
the topmost point of the fit to a parabolic shape (regions I and II in Fig. ^, respectively). 
We were motivated by Ref. to use this distance as a useful characterization since it can 
be measured in the laboratories 0. 

To perform such a fit to compute A, a problem is that transition between the behaviors 
of regions I and II is rather smooth, as shown in Fig. |[ The parabolic profile extends down 
from region I but is cut off because of the appearance of side-branches. To exclude any 
ambiguity and thus to reduce error, we have found that the quantity ^■„j(x — Xtip)/{x — Xup), 
where is the interfacial position, provides a valuable criterion to distinguish the sharp- 
cornered tip from the parabola. In region I, it is essentially hyperbolic; on the other hand, 
in region II, it reaches a constant value equal to the curvature at the tip of the parabola. 
Using this the interface belonging to region II could be extracted, and the fit done in the 
usual way. 

In Fig. 1^ the curvature as a function of the steady-state velocity is shown. We have 
obtained the same scaling behavior as that for the purely parabolic dendrites p ~ v~^^'^. 
In our case the radius p is computed using the fitted parabola discussed in the last paragraph. 
Furthermore, we have found that the values obtained for p'^v lie within 7% of the curve 
corresponding to the prediction of microscopic solvability theory [Q for 7 = 0.1. In this 
sense the tip structure, although showing clear sharp corners, acts effectively as a parabolic 
one which fits the traihng region of the dendrite. This fact suggests that the relevant quantity 
to be taken into account for the sharp-cornered tip shape is p, rather than the real curvature 
at the sharp tip, which is 10 to 20 times larger than 1/p and could be strongly conditioned 
by the lattice discretization. 

The fact that the aforementioned scaling holds, at least for the range of undercooling 
values considered here, gives a first indication that there is still only one relevant macroscopic 
length in this growth regime. A further explicit check can be done. In Fig. ^ the measured 
length A, as defined above, as well as the fitted radius of curvature p for different values of 
undercooling is shown. As A decreases (for data points on the right of the figure), the sharp 
corner at the dendrite tip is easily distinguishable and a scaling A ~ p is obtained. This 
scaling is in agreement with previously reported experimental results ||^. The whole tip 
structure is therefore self-similar including the tip region of the dendrite. For large values 
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of undercooling, the fitted p is comparable to the radius of curvature at the tip; thus the 
measure of A is no longer reliable, shown by the saturation of its values on the left part of 
Fig. |. 

The appearance of sharp-cornered tips as the steady-state velocity is reduced is not 
abrupt but gradual. This suggests that the transformation from a dendrite with a parabolic 
tip profile and a strong side-branching process in the higher velocity regime, to a sharp- 
cornered tip structure as well as "faceted" looking side-branches whose growth is rapidly 
inhibited (Fig. |I]), is a crossover behavior and not a dynamical transition. This result could 
be interpreted in terms of a kinetic roughening phenomenon [^,^ where the pattern with 
sharp-cornered tips at low driving force is kinetically roughened at large driving force, leading 
to the usual parabolic tips. To explore this idea, we have looked the change in morphology 
of the dendrite depending on the surface tension anisotropy 7, at a fixed undercooling value. 
The inset in Fig. | shows Qup which is the corner angle of the dendrite tip, or more precisely 
the extrapolation of the two sides around the tip region to function of 7. We 

observe a gradual increase of Qup asymptotically to 45° which is the angle that, according 
to the Frank diagram [|l^] and the Wulff construction , minimizes the equilibrium surface 
free-energy when 7 = 1. This energetic consideration seems play a role at the high anisotropy, 
or similarly, lower undercooling regimes. In this regime the nonequilibrium characteristics 
of the evolution appear to become less important relative to the equilibrium requirement of 
free energy minimization. 

In summary, we have studied the morphology appearing in dendritic growth of a su- 
percooled pure liquid in the low driving force regime. We found that dendrites develop 
sharp-cornered tips and the growth and coarsening of the side-branches are inhibited at 
very low undercooling. The scaling behavior for growth at large velocities, p^v ~constant, 
is still valid at low undercoolings if we use the fitted p as described above. Furthermore, 
this value of p'^v is quantitatively consistent with the microscopic solvability theory. We 
found that there is only one relevant macroscopic length scale in this slow growth regime, 
such as the fitted p. The length associated with the size of the sharp corner scales with this 
fitted p as A ~ p, and this is consistent with experimental measurements. The appearance 
of the sharp corner is gradual as undercooling decreases or anisotropy increases. Finally 
we point out that our work focuses on the sharp-cornered tip region of the dendrite, and 
the steady-state shapes, whereas in the experiment of Ref. true facets were apparent, 
suggested that it was carried out below the roughening transition of the material. It is 
not clear how to treat equilibrium roughening transition within a phase-field model where 
interfaces are diffuse (with thickness ~ e), but this is an interesting problem and should be 
pursued further in order to quantitatively compare with experiments in the very slow 
growth regime below the roughening transition. Other interesting directions include detailed 
calculations of the dynamics of dendritic growth in this low undercooling regime, and the 
investigation of this regime in three spatial dimensions. We hope to report on these studies 
in the near future. 
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FIGURES 

FIG. 1. Interfacial profile for A = 0.25 at a time t = 1.00. The growth of side-branches is 
inhibited rapidly after their appearance. 

FIG. 2. The sharp-cornered tip region of a dendrite, for A = 0.07 at a time t = 3.50 (squares). 
Also shown is the definition of A, the corresponding parabolic tip (continuous line), as well as the 
corresponding regions I and II. 

FIG. 3. The radius of curvature p of the associated parabola of the dendrite as a function of 
the velocity for different values of A (circles). The continuous line corresponds to a p ~ v^^^'^ fit, 
while the dashed line is the prediction of microscopic solvability theory for 7 = 0.1. The error bars 
correspond to the numerical error in the parabolic fit of the interfacial profile. 

FIG. 4. The distance A of the sharp-cornered tip to the fitted parabola as a function of 
the radius of curvature. The linear fit (continuous line) is extrapolated to the origin (dashed line). 
Inset: the corner angle Qup of the dendrite tip as a function of the surface tension anisotropy 
parameter 7, for a fixed undercooling A = 0.10 
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